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ABSTRACT 

What is the relevance of major mergers and interactions as triggering mechanisms for active galactic 
nuclei (AGN) activity? To answer this long-standing question, we analyze 140 XMM-Newton-se\ected 
AGN host galaxies and a matched control sample of 1264 inactive galaxies over z ~ 0.3-1.0 and 
M* < 1O 117 M with high-resolution HST / ACS imaging from the COSMOS field. The visual analysis 
of their morphologies by 10 independent human classifiers yields a measure of the fraction of distorted 
morphologies in the AGN and control samples, i.e., quantifying the signature of recent mergers which 
might potentially be responsible for fueling/triggering the AGN. We find that (1) the vast majority 
(>85%) of the AGN host galaxies do not show strong distortions, and (2) there is no significant 
difference in the distortion fractions between active and inactive galaxies. Our findings provide the 
best direct evidence that, since z ~ 1, the bulk of black hole (BH) accretion has not been triggered 
by major galaxy mergers, therefore arguing that the alternative mechanisms, i.e., internal secular 
processes and minor interactions, are the leading triggers for the episodes of major BH growth. We 
also exclude an alternative interpretation of our results: a substantial time lag between merging and 
the observability of the AGN phase could wash out the most significant merging signatures, explaining 
the lack of enhancement of strong distortions on the AGN hosts. Wc show that this alternative scenario 
is unlikely due to: (1) recent major mergers being ruled out for the majority of sources due to the 
high fraction of disk-hosted AGN, (2) the lack of a significant X-ray signal in merging inactive galaxies 
as a signature of a potential buried AGN, and (3) the low levels of soft X-ray obscuration for AGN 
hosted by interacting galaxies, in contrast to model predictions. 

Subject headings: galaxies: active — galaxies: evolution — galaxies: interactions — quasars: general 
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There is a general agreement that supermassive black 
holes (BHs) lie at the centers of nearly all galaxies, or 
at least those with a bulge component. Additionally, 
strong correlations exist between the BH mass and 
various properties of the galactic bulg e , including lumi- 
nosit y (jKormendy fe Richstond 119951: iMagorrian et al.1 
19981). stellar veloci t y dispersion (iGebhardt et al.l 12000b 
Ferrarese fc Merrittl T2000b iTremaine et all I2002D. and 



stellar mass (|Marconi fc Hunt] 120031: lHaring fcRixl 
120041) . While it has been recently proposed that 
these correlations are just the product of a statistical 
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conve rgence of s e veral galaxy mergers over cosmic 
time (penal l2007t Uahnke fe Macciol l20Tol) . these cor- 
relations have often been interpreted as the signature 
of coupl ed evolution between t he BH and its host 
galax y jKauffmann & HachneltJ |200Ct iVolonteri et al 
2001 iWvithe fe Loebl I2003L [Granato et al.1 1200 



Hopki ns et al.ll2007t iSomerville et al.ll2008D . 



Given that most galaxies are believed to have under- 
gone a quasar phase, and that the central BH represents 
a reli c of this event (jLvnden-Belll Il967t IRichstone et al.l 
1998), the co-evolution picture is naturally very appeal- 
ing even while some aspects of it remain unclear. It 
has been suggested that most of the mass of the BH is 
built up during the brightest periods of this quasar phase 
(|Soltanl [19821 lYu fe Tremainel [20021) . If there is such a 
connection between the growth of the BH and its host 
galaxy, periods of quasar activity should occur alongside 
the growth of the bulge, and the mechanism that trig- 
gers the accretion onto a once quiescent BH, turning it 
into an active galactic nucleus (AGN), should be tightly 
linked with the overall evolution of the host galaxy. The 
nature of AGN triggering is therefore of key importance 
for our understanding of galaxy evolution in general. 

According to the current paradigm of hierarchical 
structure formation, major mergers are a crucial element 
in the assembly and growth of present-day galaxies (e.g 



Kauffmann et al.l[l99a i Cole et a l. 2000; Somerville et al 
2001b IBell et al.M2006t Ijogee et al.ll2009t IRobaina et al 
2010). A closer look into the behavior of simulated col- 



lisions b etween galaxies, beginnin g with the pioneering 
work of iToomre fe Toomrel ([1972D . suggests that grav- 
itational interactions are an efficient way of transport- 
ing material toward the very center of a galaxy. Merg- 
ers and strong interactions can induce substantial grav- 
itational torques on the gas content of a galaxy, de- 
priving it of its angular momentum, leading to inflows 
and the buildup of huge reservoirs of gas in the cen- 
tcr (iHernquistl [19891: IB arncs fe Hernauisd U991 U9961; 
Mihos fe HernauistJll996llSpringel et al II2005I; I Cox et al.l 
20061 iDi Matteo et alJl2007t ICox et alJl2008[ ). 

From early on, major mergers have been related to ob - 
servations of powerful nuclear starbursts (jGunnl IT979I ). 
and c onnections with quasar activity were made soon 
after. iStocktonl (H982), in a study of luminous quasars 
with close companions, suggested that these neighbor- 
ing galaxies could be survivors of a strong interaction 
with the quasar. Further observational studies came to 
support this picture: more cases of quasars with close 
companions were found, and post-merger features were 
detected in the host galaxies, whenever it was possible 
to res olve them (e.g..lHeckman et al.lll984HGehren et al.1 
19841: [Hutchines et al.lll984L [19881 : iStockton fe Ridgwavl 
19911 : iHutchings fe Ncff 1992). The merger-quasar con- 



nection scenario gained strength with the discovery of the 
ultraluminous infrared galaxies (ULIRGs). More than 
95% of these were found in a merging state, some of 
them hosting an AGN. This led to the scenario in which 
ULIRG s and quasars were p art of the same chain of 
events dSanders et al.lll988alfbl: ISanders fe Mirabellll996l: 
Surace et al.lll998t iSurace fc SandersHl999T Ts urace et aLl 
2000; [Canalizo fe Stocktonll2000l |2001[ ). 



Many observational studies of luminous AGN found a 
high rate of merging signatures in their hosts and de- 
tected the presence of very close comp anions, which 
before HST could not be resolv ed Ce.g.. iBahcall et all 
19971 iCanalizo fe Stockton! 120011: iZakamska et al.1 120061 : 
Urrutia et al.1 [20081 ) . At the same time, deeper imag- 
ing of AGN host galaxies that were initially classified 
as undisturbed revealed post-merger feat ures not previ- 
ously detected, both from space-based ( Canaliz o et al.l 



ously 
2007; 



With the advent of the Hubble Space Telescope (HST), 
deep imaging of AGN host galaxies at higher red- 
shifts became possible with unprecedented resolution. 



Bennert et al.l l2008|) and g round-based observa- 
tions ( Ramos Alm eida et aLll2010l ). 

There is, however, one major caveat for most of the 
studies listed above: almost none of them made use of, 
or had the access to, an appropriate control sample of in- 
active galaxies; such a control sample is essential for dis- 
cerning if the merger rate is in fact enhanced with respect 
to the "backgr ound level" , i.e . , the m erger rate of inactive 
galaxies. Only iDunlop et al.l (|2003| ) compared their sta- 
tistically complete sample of quasars against the quies- 
cent galaxy population, finding no difference in the struc- 
tural parameters between samples, as well as no enhance- 
ment in the large-scale disturbances. Even if not explic- 
itly, this showed a clear divergence from previous stud- 
ies regarding the merger- AGN connection scenario, and 
agreed with the very low frequency of post-merger sig- 
nature s observed on Seyfert galaxies and low- luminosity 
AGN (IMalkan et al.lfl998t ISchade et al.ll2000T) . 

A new era of large HST programs now offers the po- 
tential for resolving this discrepancy. The imaging of 
larger, contiguous fields has yielded a large number of 
objects, making it possible to study AGN hosts at space- 
based resolution, and at the same time to compile a con- 
trol sample of no n-active galaxies. In itial studies using 
HST imaging by iSanchez et all (|2004l ) with the Galaxy 
Evolution from Morpholo gies and SED s surve y (GEMS, 
iRix et al.l I2004T ) and by fcrogin et all ()2005l ) with the 
Great Observatories Origins Deep Survey (GOODS, 
Giavalisco ct al. 2004) found no evidence for an enhance- 
ment in merging signatures of AGN hosts over control 
galaxy samples. If merger activity does not play a major 
role in AGN triggering, other methods to produce gas in- 
flows, build up the bulge, and fuel the BH should also be 
of importance. Alternate secular mechanisms — minor in- 
teractions, large scale bars, nuclear bars, colliding clouds, 
supernova explosions — can also lead to angular momen- 
tum removal and gas inflows fr om different scales to the 
central regions (for reviews, see lKormendy fe Kennicuttl 
l200i I Wadal2004t IMartinl200i IJogeeEjo 6J) . While these 
processes have usually been relate d to Seyfert galax- 
ies and low-luminosity AGN (e.g.. iSimkin et al.lll980l : 
iTaniguchil [19991 Hopkins fe Hernquistll2009D . thev could 
potentially play a larger role than usually reckoned for 
more luminous AGN as well. Although the results from 
the GEMS and GOODS surveys are highly intriguing, 
the field sizes of ~0.22 deg 2 and ~0.08 deg 2 respectively 
were still too small for definitive conclusions to be drawn. 
A suitably larger sample would be required to turn these 
appealing hints into statements. 

In this context, we tackle this long-standing issue 
by performing a comprehensive morphological analysis 
of a sample of X-ray-selected AGN host galaxies from 
the Co smic Evolution Survey (COSMOS, iScoville et aLl 
12007b ) , the largest contiguous area ever imaged wi th the 
HST (jScoville et al.|[2007al; iKoekemoer et alj |2007). Our 
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goal is to disentangle the actual relevance and predomi- 
nance of major galaxy mergers from the other suggested 
mechanisms for the fueling of the BH. 

In the past, targeted high- resolution imaging of AGN 
hosts has only been possible for small samples, while ex- 
tensive ground-based surveys with large samples have 
lacked of the necessary resolution to perform detailed 
morphological studies at moderate redshifts. Ea r lier re - 
sults from the detailed analysis by iGabor et al.l (2009), 
where the morphologies of ^400 AGN host galaxy can- 
didates from the COSMOS field were parameterized, 
showed that these had an asymmetry distribution con- 
sistent with that of a control sample of inactive galax- 
ies, and lacked an excess of companions, already sug- 
gesting that major interactions were not predominant 
among AGN as a triggering mechanism. Here we use 
the largest sample of optically confirmed X-ray-selected 
AGN ever imaged at HST resolution from the COSMOS 
survey and perform a visual inspection of the morpholo- 
gies of the host galaxies. We opt for a visual analysis 
of our galaxies over an automatic classification system 
because of the inherent problems and incompleteness of 
the latter in identifying mergers, even for some obvious 
case s, as cautioned by recent studies probing both meth- 
ods (|Jogee et al.ll2009HKartaltepe et alj|2010j) . To estab- 
lish the relevance of our findings, we compare the AGN 
hosts to a matching sample of inactive galaxies from the 
same exact data set. 

Throughout this paper we assume a flat cosmology 
with H = 70 km s" 1 Mpc -1 , Cl M = 0.3, and n A = 
0.7. All magnitudes are given in the AB system unless 
otherwise stated. 

2. DATA SET AND SAMPLE 

We will perform our analysis on a complete sample 
of X-ray selected optically confirmed type-1 and type-2 
AGN from the COSMOS field. 

The COSMOS survey features the largest contiguous 
area ever imaged with the HST. The location of the 1.64 
deg 2 field, close to the celestial equator, allows access 
from several major space and ground-based observato- 
ries, enabling a large multiwavelength coverage from X- 
ray to radio from supplementary observational projects. 

One of the most effective ways of finding AGN is to 
make use of the X-ray em ission due to the accreting 
BH (e.g.. lMushotzkvll200l . Complete coverage of the 
whole COSMOS field in X-rays was achieved wi t h the 
XMM-Newton (XMM -COSMOS, IHasinger et all 120071: 
ICappelluti et al Il2009l ) through 55 pointings with a total 
ex posure time of ~L5 Ms. We use the catalog presented 
in IBrusa et al.l (|2010l ). which provides the most likely 
optical and infrared counterparts to the X MM sources 
based on a likelihood ratio technique (see IBrusa et al.l 
2007 for details). 

From the X-ray catalog we draw a parent sample 
of ^550 sources c lassified as type-1 AGN from spec- 
trosc opic surveys (|Trump et al.l 120071 120091: iLillv et al.l 
120071) revealing broad emission line s, and from spec- 
tral energy distribution (SEP) fitting (jCapak et al.1120071: 
ISalvato et aLll2009l: [iTbert et al.i r20091. We also include a 
subsam ple of X-ray s e lected type-2 AGN, based on those 
used by IGabor et al.l ([2009) drawn fr om a parent sample 
of ~3 00 narrow emission line objects ([Trump et al.ll200~7l 
[20091) . 




43 44 
log L x (2-10 keV) (erg/s) 

Figure 1. The X-ray luminosity distribution of our sample in the 
2-10 keV energy band (solid line). For reference, we also show the 
distribution of the type-1 AGN subsample only (dotted line). 



In this paper, we analyze the morphological proper- 
ties of the AGN host galaxies. For this, we take ad- 
vantage of the high-resolution imaging of the COSMOS 
field with the HST. These observations comprise 583 
orbits using the Advanced Camera f or Surveys (ACS) 
with the F814W (broad /-band) filter (jKoekemoer et al] 
120071 ) . The imaging data feature an oversampled scale of 
O'.'03/pixel. Although the ACS survey of the COSMOS 
field is highly homogeneous, the exact depth achieved is 
dependent on the angle of the telescope with the Su n 
at the time of the observations ([Leauthaud et al1l2007|) . 
Ninety six out of the 575 pointings were made with an 
angle smaller than a critical value of 70°, leading to a 
slightly shallower image. The limiting surface brightness 
levels above the background for the pointings made with 
an angle with the Sun larger and smaller than the critical 
value are ~23.3 mag arcsec -2 and ^22.9 mag arcsec -2 
respectively. 

We restrict our sample to the rcdshift range z ~ 
0.3 — 1.0. For the majority of our final sample, we used 
high-confidence spectroscopic redshifts, wh ile for the rest 
(20%) , we used photometric redshifts by ISalvato et al.l 
(2009). The lower redshift cut is chosen due to the low 
number of AGN below z ~ 0.3 (see Section 3 for fur- 
ther details), and also to avoid working with saturated 
sources. The upper limit arises because the F814W filter 
is shifted into rest- frame UV for sources above z ~ 1. 
This would mean that we would be specifically looking 
at the light from young stars and star formation knots, 
and therefore at biased morphologies. At the same time, 
for the case of the type-1 AGN, the bright nucleus would 
start to dominate due to its blue color, strongly outshin- 
ing the host and making it nearly impossible to resolve. 

Given our interest in the morphologies of our sample of 
AGN host galaxies, we decided not to consider galaxies 
fainter than Ifsuw — 24. Visual morphological classi- 
fication of these objects would be particularly difficult, 
and we determined that no consistent information could 
be extracted at this magnitude. For the case of the type- 
1 subsample, we applied this criterion after the nucleus 
removal (see Appendix A for details). 

This selection yields 83 type-1 and 57 type-2 AGN. The 
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median redshifts of the type-1 and type-2 subsamples are 
0.80 and 0.67 respectively. The type-2 subsample has a 
median apparent Ifsuw of 20.9, slightly brighter than 
the type-1 subsample with 21.7 (after nucleus removal, 
see Appendix A) due to the higher median redshift of the 
latter. 

Figure Q] shows the X-ray luminosity distribution of 
our sources in the 2-10 keV energy band. The val- 
ues were obtaine d mainly from those calculated by 
Lusso et al.l (1 20 1 Of) an d are complemented with those by 
Mainieri et al.l ()2007l ). The median of our distribution 
lies at Lx = 10 43 5 erg s _1 , which means that we are 
probing a reasonably luminous representative AGN sam- 
ple. For reference, in Figure [1] we also show the X-ray 
luminosity distribution of the type-1 AGN subsample 
only, which dominates the overall distribution and has 
a slightly brighter median Lx (10 43 6 erg s _1 ) than the 
type-2 subsample (10 43 3 erg s" 1 ). 

3. METHODOLOGY 

In this paper we analyze the morphologies of a sample 
of AGN host galaxies and of a control sample of inactive 
galaxies using high-resolution HST/ ACS single /-band 
images. In the following subsections we explain how we 
built the comparison sample (hereafter CS), the motiva- 
tion of choosing a visual inspection over an automatic 
method, and the classification scheme used. 

Analyzing the host galaxies of type-1 AGN is complex, 
due to the presence of the bright active nucleus in the 
images, that, depending on the contrast, can outshine the 
host galaxy to different extents. We overcome this issue 
through a two-dimensional decomposition of the AGN 
and its host galaxy, modeling the bright active nucleus 
with a point-spread function (PSF) and the host galaxy 
with a iSersid (1968) profile using GALFIT (|Peng et all 
[20021 120101) . and then removing the nuclear contribution. 
This process is described in detail in Appendix A. 

3.1. Comparison Sample 

The large number of galaxies available from the COS- 
MOS HST observations provides us with the unique op- 
portunity of building a control sample from the same 
data set that we draw our AGN from. For our study, 
we require the control sample to permit us elaborate a 
comparison regarding distortion features. On this re- 
spect the most relevant parameter is the signal-to-noise 
ratio (S/N), hence we construct the comparison sample 
by selecting inactive galaxies matching each AGN both 
in apparent Ifsiaw magnitude and photometric redshift. 
This is both required and sufficient since (1) the S/N de- 
termines the visibility of the merger signatures, and (2) 
while the stellar masses might differ slightly (factors of 
~2 at the same distance and brightness), the mass de- 
pendence on the merger rate is no t strong, with only a 
modest increase for higher masses ()Bundv et al.l l2009). 

Specifically, for each AGN host galaxy we select 10 sim- 
ilar comparison galaxies fro m the COSMOS ACS cata- 
log (|Leauthaud et all 12007). Each selected comparison 
galaxy is required to have an Ifsuw magnitude within 
a range of AIpsuw =0.1, and a photometric redshift 
within a range of Az = 0.05. If not enough galaxies were 
found, the search ranges were increased by 10%. On av- 
erage, 1.8 iterations were performed to find the required 



number of inactive galaxies for each AGN host. For the 
case of the type-1 AGN subsample, the magnitudes of 
the host galaxies after the removal of the active nucleus 
are used for the selection of the control sample. 

With the inactive comparison sample in hand, we re- 
move galaxies that are unlikely to be AGN host galaxy 
counterparts a priori via an initial visual inspection. 
Such galaxies include: (1) bulgeless disks and irregulars, 
which would represent a low-mass population, having no 
corresponding partners on the AGN sample, and (2) for 
the type-1 subsample, edge-on disks, which could in prin- 
ciple hold an AGN but this would be heavily obscured 
and therefore not be a type-1. 

Finally, the construction of the control sample for the 
type-1 AGN requires an additional effort. The nucleus re- 
moval process usually leaves residuals in the center which 
certainly affect any blind classification, making the type- 
1 AGN host galaxies readily discernible from the control 
sample. To resolve this problem, we mock up our selected 
inactive galaxies as AGN by adding a star in the center 
as a fake active nucleus, as we describe in Appendix B. 
We then apply the same subtraction procedure as for the 
original type-1 AGN, attempting to make the two sam- 
ples indistinguishable. As we show in Appendix C, any 
effects on the selection of the comparison sample due to 
flux variations caused by the nucleus subtraction process 
can be neglected. 

Our final comparison sample consists of 1264 galaxies 
in total. The Ifsuw an d redshift distributions of the re- 
sulting type-1 and type-2 comparison samples are consis- 
tent with being drawn from the same parent distribution 
as the AGN subsamples, even after the removal of the 
unlikely counterparts described above. A Kolmogorov- 
Smirnov test on each couple of Ifsuw and redshift distri- 
butions confirms with probabilities > 38% that the AGN 
and control samples are consistent among each other (in 
general < 5% is used to show that two distributions dif- 
fer). 

3.2. Visual Classification 

Merger events come in many different flavors due to the 
large parameter space involved (e.g., merger stage, view- 
ing angles, mass ratio, and gas fractions). Sometimes 
they can be obvious at first sight, but some others can 
be very subtle, or simply undetectable at the sensitiv- 
ity of the observations. At our redshift range and image 
resolution, it has been shown that automatic classifica- 
tion methods to identify mergers tend to miss several 
obvio us cases, and cannot compete with visu al inspec- 
tion (jJogee et al.l 120091 : iKartaltepe et alll2010D . On the 
other hand, when the numbers involved are over the tens 
of thousands, visual classification becomes impractical 
and an automatic approach would be needed. General 
measurements of structural parameters that can be cor- 
related with some physi cal process have prov en to be a 
good compromise (e.g.. iReichard et "all 12009. using the 
lopsidedness as a tracer of merging and star formation). 

Considering the above, in this paper we opt to identify 
merger and interaction signatures visually. The number 
of objects we are dealing with allows us to do so (~1400 
in total), and the image quality deserves a detailed case- 

3 Wit h the notable exce ption of the citiz en-based G alaxy Zoo 
project (Lintott ct al. 2008, http://www.galaxyzoo.org). 
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Table 1 

Results from the Visual Analysis by the 10 Classifiers. 



Classifier 


MC 


KI 


KJ 


JK 


AK 


TL 


AR 


MS 


KS 


JT 




^ AGN 


140 


140 


140 


140 


22 


40 


57 


140 


38 


98 




Ntype-1 


83 


83 


83 


83 


22 


10 





83 


19 


11 




Ntype-2 


57 


57 


57 


57 








57 


57 


19 


57 




Ncs 


1264 


1264 


1264 


1264 


177 


357 


537 


1264 


357 


903 


- 


Hubble type 
























Bulge AGN 


25.7% 


51.4% 


31.4% 


20.0% 


40.9% 


55.0% 


29.8% 


43.6% 


26.3% 


37.8% 


35.2% ± 11.0% 


DiskAGN 


74.3% 


48.6% 


68.6% 


80.0% 


59.1% 


45.0% 


70.2% 


56.4% 


73.7% 


62.2% 


64.8% ± 11.0% 


Bulge cs 


24.6% 


43.3% 


29.6% 


25.2% 


47.5% 


54.3% 


29.4% 


43.6% 


17.9% 


40.5% 


34.370 ± 9.5%> 


Disk cs 


75.4% 


56.7% 


70.4% 


74.8% 


52.5% 


45.7% 


70.6% 


56.4% 


82.1% 


59.5% 


65.7% ± 9.5% 


Distortions 
























Dist-O^GAT 


62.9% 


43.6% 


48.6% 


56.4% 


50.0% 


47.5% 


71.9% 


56.4% 


47.4% 


55.1% 


54.2% ± 7.5% 


Dist-Ocs 


65.5% 


47.3% 


60.1% 


63.0% 


67.8% 


51.5% 


78.0% 


58.5% 


51.5% 


61.2% 


59.9% ± 7.6% 


A-Dist-0 


-2.6% 


-3.7% 


-11.6% 


-6.5% 


-17.8% 


-4.0% 


-6.1% 


-2.1% 


-4.2% 


-6.1% 


-5.6% ± 3.5% 


Dist-lAGiV 


24.3% 


26.4% 


45.0% 


32.9% 


40.9% 


40.0% 


21.1% 


30.0% 


50.0% 


16.3% 


30.8% ± 9.3% 


Dist-l C s 


22.4% 


28.3% 


34.2% 


26.9% 


19.2% 


34.2% 


16.0% 


33.5% 


39.2% 


17.8% 


27.5% ± 6.5% 


&Dist-l 


1.9% 


-1.9% 


10.8% 


6.0% 


21.7% 


5.8% 


5.0% 


-3.5% 


10.8% 


-1.5% 


3.2% ± 5.7% 


Dist-2^ Gi v 


12.9% 


30.0% 


6.4% 


10.7% 


9.1% 


12.5% 


7.0% 


13.6% 


2.6% 


28.6% 


15.0% ± 8.8% 


Dist-2 C s 


12.1% 


24.4% 


5.7% 


10.1% 


13.0% 


14.3% 


6.0% 


7.9% 


9.2% 


20.9% 


12.6% ± 6.5% 




0.8% 


5.6% 


0.7% 


0.6% 


-3.9% 


-1.8% 


1.1% 


5.7% 


-6.6% 


7.6% 


2.4% ± 3.5% 



Note. — Wc indicate the number of objects classified by each person for the AGN sample and individually for each type-1 and typc-2 subsamplc, 
as well as for the comparison sample (CS). For the distortion classifications, we include the difference between samples as A^ist-x— Dist— Xagn — 
Dist— Xc"S- For each category, wc include the mean, and its dispersion, weighted according to the number of objects classified by each person. 



Table 2 

Comparison of our mean Hubble type classification 
with that of parametric estimators of galaxy 
morphologies 







GALFIT b 


ZEST C 


Bulge^GiV 


35.2% 


25.7% 




Disk^GiV 


64.8% 


55.0% 




Bulgecs 


34.3% 


41.2% 


19.6% 


Disk C s 


65.7% 


43.5% 


67.8% 



a Weighted mean of the 10 classifications (as in Table FIT), 
k Percentages over 100% of the samples. The rest of the 
objects had an intermediate Scrsic index (with 2 < n < 3, 
sec Appendix A for details). 

c Percentages over 100% of the comparison sample. The 
remaining galaxies were classified as irregulars (11.3%) by 
ZEST, and a few did not make it to the catalog (1.2%). 

by-case examination. 

These visual studies can be subjective. In our case, the 
absolute fraction of merging galaxies measured by visual 
classifiers will depend on their own experience and back- 
ground, and hence it is plausible that they can differ 
substantially among each other. Nevertheless, any per- 
sonal scale and criteria each classifier uses will be applied 
equally on both samples, active and inactive galaxies. 
Therefore, a key quantity on our study will be, more than 
the absolute fractions of merging galaxies, the difference 
between the merging fractions measured by a given classi- 
fier. If we instead focus on how each individual classifier 
perceives one sample compared to the other, by consid- 
ering the differential between the merging fractions of 
active and inactive galaxies, this subject iveness can be 
accounted for. Furthermore, the consistency of this study 
is improved by (1) using ten independent human classi- 
fiers to add statistical robustness and (2) mixing both 
samples of active and inactive galaxies so that the clas- 
sification is actually blind and therefore does not favor 
either the AGN hosts or the inactive galaxies. 

We break the classification down into two parameters. 

I. Hubble type. We attempt to state whether the host 



galaxy belongs to one of the two basic morpholog- 
ical classes: bulge or disk dominated. 

2. Distortion class. We define three classes regarding 
the degree of distortion of the galaxy as follows. 

(a) Dist-0. Galaxies that appear undisturbed, 
smooth and/or symmetric, showing no inter- 
action signatures. This also applies to cases 
where the small diameter of the galaxy does 
not allow a detailed analysis. We pay par- 
ticular attention to self-induced asymmetries 
such as dust lanes or star-forming regions, 
which are usually seen as small clumps in well- 
resolved spirals. 

(b) Dist-1. Here we include objects with mild dis- 
tortions. This could be due to a minor merger 
for example, but at the same time could also 
be because of low S/N. This interaction class 
is a "gray zone" in which most of the discrep- 
ancies in the classification between the 10 peo- 
ple arise. 

(c) Dist-2. Strong distortions, potential signs for 
ongoing or recent mergers. This class mainly 
includes galaxies which have highly disturbed 
morphologies or show visible signatures of 
strong interactions, such as large tidal tails, 
arcs, debris, etc. Double-nucleus systems also 
fall into this category. 

Illustrative examples of the Hubble type and distortion 
classes are shown in Figures [2] and |3l respectively. 

For the visual inspection, the classifiers had access to 
FITS images which they could re-scale in order to look 
for high-contrast and subtle features that may have not 
showed up at an arbitrary brightness scale. 

4. RESULTS 

The results from the visual classification by 10 people 
(MC, KI, KJ, JK, AK, TL, AR, MS, KS, and JT), for 




AGN host galaxies 



Inactive galaxies 



< 


> 











Figure 3. Example galaxy images arranged into different distortion classes with 100% agreement between the independent classifiers. 
The cutouts are 4'.'8 X 4'.'8. Black residuals at the center of some of the galaxies are residuals from the point source removal. 



both Hubble type and distortion classes, are shown in 
Table [TJ For the different distortion classes, we show the 
difference between samples (hereafter A) as the distor- 
tion fraction of the AGN minus that of the control sam- 
ple. The results are weighted according to the number 
of objects classified by each persorQ and used to calcu- 
late the mean fractions, /i, which we also display in Ta- 
ble [TJ Figures [2] and [3] show examples of active and inac- 
tive galaxies which were classified with 100% agreement, 
arranged into the different Hubble type and distortion 
classes respectively. 

4 Each classifier looked at a minimum of ~200 galaxies from the 
combined sample. For each classifier the samples were shuffled, to 
assure that even if all of them decided to look at 200 galaxies, they 
would be looking at different objects. On average, each galaxy was 
classified 6.3 ± 1.0 times. 



4.1. Perception of the Hubble type 

No morphology priors are applied in the selection of 
our comparison sample, with the minor exception of the 
pruning of irregulars and edge-on disks as described in 
Section 3.1. In order to test whether the samples are 
consistent regarding their morphological composition, we 
compare the AGN and comparison samples in Table [TJ 
Although the mean values show a high dispersion due to 
the large discrepancies between classifiers, the results for 
the AGN and comparison samples are in good general 
agreement for each classifier. 

The large fraction of disks, in particular in the AGN 
sample, is interesting. To verify if this could be due to 
systematic bias by the classifiers, we will use two inde- 
pendent parametric estimators of the morphological type 
available at hand. First we compare the results from our 
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Figure 4. Distributions of the difference in the Monte Carlo sampled distributions of Dist-2 fractions between the AGN and control 
samples for the ten classifiers. For each distribution, a deviation from zero difference (dotted line) towards positive values indicates a higher 
fraction of distorted active galaxies, whereas a deviation towards negative values shows a higher fraction of distorted inactive galaxies. 



GALFIT models chosen earlier, which we extended to 
our type-2 subsample as well as to its comparison galax- 
ies. We identify sources as bulge- or disk-dominated if 
the best-fit results from GALFIT had Sersic indices of 
n = 4 and n = 1 respectively. The rest of the galaxies 
fell between the two. As a second test, we look up the 
results for our comparison sample from the Zurich Es- 
timat or of Structural Types (ZEST, see iScarlata et al.l 
120071 for details), in which the structure of thousands 
of COSMOS galaxies was quantified through a principal 
component analysis over a combination of Sersic index 
and five non-parametric diagnostics. The ZEST results 
show the fractions of galaxies classified either as bulges or 
disks. Of the remaining fraction classified as neither, the 
majority (11.3%) was classified as irregulars, most likely 
due to the lack of sensitivity of these automatic classi- 
fication schemes to peculiar systems such as interacting 
galaxies; this is consistent with the observed fraction of 
highly distorted comparison galaxies. Sixteen galaxies, 
accounting for the remaining 1.2%, did not make it into 
the catalog. 

Table [2] shows both of these tests along with the 
weighted mean fractions for comparison. It is clear that 
the numbers from these tests follow the trend seen in the 
visual classification. These tests provide a lower limit to 
the fraction of disks, with > 55% of our AGN sample 
being hosted by true disks. 

4.2. The distortion fractions 

Our prime interest lies in the observed difference in dis- 
tortion fractions between samples of active and inactive 
galaxies. The absolute values in distortion fractions de- 



termined by the 10 classifiers are of lesser interest since 
the internal calibration for the three distortion classes 
differs between the classifying individuals. Since any 
subjectiveness will be applied equally to both active and 
inactive samples, using the differences in the distortion 
fractions instead of absolute levels removes the person- 
to-person calibration differences and allows an unbiased 
interpretation. 

Considering that the merging signatures we were look- 
ing for could sometimes be faint and weak, we address 
the potential loss of sensitivity to such features due to 
the slightly shallower limiting magnitudes for ~17% of 
the pointings (i.e., those with Sun-angles of <70°). For 
each person, we have also carefully analyzed the results 
by dividing their classified sample into sources with sun- 
angles either side of this critical angle. We find that there 
is no statistically significant difference in the distortion 
fractions as a function of Sun-angle. In addition, as the 
assignment of individual objects to either a deep or shal- 
low field is effectively random, and given that the AGN 
distortion fractions are compared directly with those of 
a comparison sample selected from the same data set 
(and thus with the same limiting surface brightness is- 
sues) , the overall impact on our results of any bias toward 
smaller distortion fractions in the shallower fields would 
in any case be negligible. 

The objects that fell into the Dist-2 class were those 
which presented the strongest distortions, and hence sig- 
natures of a major interaction, for each individual clas- 
sifier As any difference in recent major merger incidence 
would show itself in this class, we will focus on the Dist-2 
results here. 
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Figure 5. Combined posterior probability distribution of the dif- 
ference of highly distorted galaxies between the AGN and control 
sample for the 10 classifiers. The central 68% confidence level is 
marked with vertical dashed lines, which shows that the histogram 
is consistent with zero difference (dotted line), ruling out any sig- 
nificant enhancement of merging signatures on our sample of AGN 
hosts with respect to the comparison sample of inactive galaxies. 



100% 



£ 60%- 




Figure 6. Cumulative distribution of the simulated Dist-2AGN 
fractions, showing the 68% and 95% confidence levels with the 
dashed lines. As mentioned in the text, this confirms with a 95% 
confidence that the highly distorted AGN fraction can not be larger 
than 24.08%. 



4.2.1. Combining 10 classifications 

In table [T] we have already listed the Dist-2 fractions 
for all classifiers, their mean values, and also the mean of 
the difference in Dist-2 fractions between the AGN and 
comparison samples. This permits the following initial 
assessment under the assumption of Gaussian errors: the 
difference (2.4%) is below the uncertainty of 3.5%, and 
hence it is not significant. 

Nevertheless, since the error distribution is in fact not 
Gaussian but follows a binomial distribution (according 
to the number of distorted galaxies in a sample of given 
size) it is important that we use the correct combination 
of results in order to give answers to the two main ques- 
tions: (1) Is there a genuine difference between the frac- 
tions of strongly distorted AGN hosts and inactive galax- 
ies, and (2) with the given sample size, what difference 
in distortion fractions between samples can we actually 
rule out at a given confidence level — in this case we chose 
95%. The first question asks whether the given dataset 



shows an enhanced AGN distortion fraction or not. The 
second question probes the discriminative power of this 
sample, and allows us to gauge the actual importance 
of a null-result in question (1), since a decreasing sam- 
ple size means an increasing uncertainty in the distortion 
fractions and hence small samples have near zero discrim- 
inative power. 

Using the correct binomial error statistics for the dis- 
tortion fractions of AGN and inactive galaxies, we com- 
pute for each classifier the probability distribution for the 
difference Ai>; st _2. This is done in the following way: 
(1) individially for each classifier, we Monte Carlo sam- 
ple their pair of Dist-2 binomial probability distributions 
for the AGN and comparison samples separately, (2) we 
compute the difference between these randomly sampled 
values, (3) we repeat this process one million times for 
each classifier which yields 10 distributions for Aui St -2, 
(4) we normalize these probability distributions by the 
Dist-2(7s values measured by each person as shown in Ta- 
ble [U in this way removing any bias applied by each clas- 
sifier's personal scale (Figure [4]), and (5) now in "differ- 
ential" space, where we are insensitive to between-person 
scatter, we combine these 10 distributions by co-adding 
their histograms, weighted by the size of the sample each 
person classified.]. 

The resulting probability distribution is shown in Fig- 
ure [5] The histogram is fully consistent with zero differ- 
ence, as indicated by the central 68% confidence interval 
denoted by the vertical dashed lines estimated by the ar- 
eas at both ends, encompassing 16% each. This confirms 
the simple analysis from above: our study shows no sig- 
nificant difference between the fractions of strong distor- 
tions of AGN and inactive galaxies. Regarding the dis- 
criminative power of our sample, in Figure [6] we show the 
cumulative distribution of the Dist-2^GAf fraction from 
Figure O The distribution shows that with 95% confi- 
dence the distortion fraction of AGN is in any case not 
larger than the inactive distortion fraction by a factor 
of 1.9, when considered relative to the mean distortion 
level found by the 10 classifiers (12.6%). Hence, the vast 
majority of AGN host galaxies at z < 1 with the given lu- 
minosities do not show signatures of having experienced 
a recent major merger. 

4.2.2. Mass dependency 

Even if there is no overall difference between the 
fractions of highly distorted AGN and inactive galax- 
ies, it is still interesting to look at the situation in 
mass-space, and investigate the possibility that an en- 
hancement of the AGN merger fraction could be hid- 
den because we consider the sample as a whole, re- 
gardless of stellar mass. Major merging is a key ele- 
ment in the assembly and evolution of massive galax- 
ies (e .g.. iBell et al J 12005 iLin et all 120081; iBundy et al.l 
l2lMlvander Wei et alll2009l ; iRobaina et al.ll2010D . and 
in order to test if the fraction of highly distorted AGN 
host galaxies is significantly enhanced at the massive end 
(higher than ~ 10 10 5 M©), we have estimated stellar 
masses for our samples of active and inactive ga laxies. 
We use the calibration from IBell &; de Jong] (|200lD based 
on the Chabrier initial mass function. By obtaining the 

5 This represents a combined Bayesian posterior probability dis- 
tribution with sample sizes as individual priors. 
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Figure 7. The combined differences in distortions of intermedi- 
ate (10 9 ' 3 < KL/Mq < 10 10 5 ; top panel) and massive (10 10 5 < 
Mj/Mq < 10 11 ' 7 ; bottom panel) galaxies are shown. In both 
cases, the central 68% confidence levels (dashed lines) are consis- 
tent with zero (dotted line). 

y-band luminosities: 

Ly = 10 -0-i(V-4.82) (1) 

and assuming a common mass-to-light ratio from the 
rest-frame (B — V) color, we derive stellar masses in solar 
units: 

M* = io-°-™+i.305(s-v) x Ly (2) 

with all magnitudes in Vega zero point. 

For the inactive galaxies and the type- 2 subsample we 
ob tain rest-frame B a nd V from the photometric catalog 
bv lllbert et~aT1 (|2010[ ). For the type-1 AGN, however, we 
cannot use that information because it includes the con- 
tribution from the luminous AGN. Therefore, we obtain 
the rest-frame y-band luminosities from the observed 
Ifsuw after the nucleus removal process and estimate 
the color term by computing the linear regression over 
the rest-frame (B — V) colors as a function of redshift for 
the type-2 AGN. This yields the relation 

(B — V)vega = 0.136 z + 0.541. (3) 

The combined differences of highly distorted galax- 
ies for two bins of stellar mass (10 9 3 10 10 5 Mq and 

10 10 - 5 lO n - 7 M ) are shown in Figure For both 

samples the ratio of galaxies occupying the massive bin 
is roughly 2:1 relative to the less massive one, hence we 
are dealing with very massive galaxies. Even if for the 
galaxies with stellar masses higher than ~ 10 10 ' 5 M Q 
there is a modest enhancement in the distortion fraction 
of the AGN hosts over the control sample (Figure bot- 
tom panel) , it is again within the 68% confidence interval, 
i.e., it is not statistically significant. Therefore, it cannot 
be considered as an empirical proof of an enhancement 
at the massive end. 

5. DISCUSSION 

From a detailed analysis of the results of our visual 
classification we showed that the fractions of heavily dis- 
torted active and inactive galaxies are consistent within 
the central 68% confidence interval and that the Dist- 
2 fraction of AGN host galaxies is less than twice that 
of the inactive galaxies at a 95% significance level, as 



shown in Figures [5] and [51 respectively. Putting these 
findings in context, provided that the duration of merger 
signatures and the visibility of the AGN phase overlap 
with each other, this indicates that there is no evidence 
that major merging plays a key role in the triggering of 
AGN activity in our sample. But what about the pos- 
sible alternative scenario in which, in spite of a causal 
connection between merging and AGN triggering, we do 
not detect an enhancement of merger signatures in the 
AGN population due to a significant time lag between 
the interaction and the start of the AGN phase? Below 
we address this possible alternative interpretation with 
some simple tests, and discuss the implications of our 
results. 

5.1. Alternative Interpretation: Time Lag Between 
Merging and the Observability of the AGN Phase 

Appealing simulations of mergers between gas-rich 
galaxies state that the peak of star formation and quasar 
activity will occur during the final stages of the interac- 
tion, close to coalescence, withi n a more relaxe d than 
distorted bulge-like remnant (jDi Matteo et al.l 120051 : 
ISpringel et al.l 120051 ) . In these models, during the first 
passage only modest starbursts are triggered and no 
major BH accretion occurs, and therefore the galax- 
ies would not be detected as AGN. Furthermore, ad 
hoc models that inclu de obscuration in galaxy mergers 
(Hopkins et al. 2005b) predict that, beginning from the 
early stages of the interaction, the AGN is "buried" for 
^90% of its lifetime by large column densities, only re- 
vealing itself toward the end of the merger. However, all 
these models work with sub-grid prescriptions of BH ac- 
cretion and fail to spatially resolve the actual accretion 
process by several orders of magnitude. 

If there is indeed a substantial time lag after merg- 
ing prior to the AGN activity becoming detectable, then 
the strong merging signatures we attempt to find could 
have already been washed out. Moreover, if AGN are ob- 
scured as the interacting galaxies coalesce, there could be 
a "contamination" population of undetected strong BH 
activity occurring within our control galaxies undergo- 
ing a major merger. Finally, a third issue related to the 
obscuration plus time lag scenario is that the observed in- 
teractions that are occurring on a fraction (^15%) of our 
AGN host galaxies should be unrelated to the detected 
BH accretion — under the assumption of a large time lag 
we would not expect to see strong merging signatures. 

5.1.1. AGN Hosted by Disks: Not a Relic from a Major 

Merger 

In the preceding text, we raised a possible alterna- 
tive explanation for our results, that most major merg- 
ers could be missed because the time lag between merg- 
ing and the observed AGN episode could be substantial, 
washing out the signatures that the HST / ACS resolution 
allows us to detect. 

Models can provide us with some clues about the ob- 
servability timescales during an inter action. For exam- 
ple, simulations of major mergers by lLotz et all (2008) 
quantified that the strong signatures could still be de- 
tected 0.7 Gyr after the merger, by degrading their snap- 
shots to the resolution of HST z ~ 1 imaging. Thus, in 
order to explain the observed zero distortion enhance- 
ment, a lag of at least 0.7 Gyr between coalescence and 
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the visible phase of the AGN would be required for all 
galaxies^. It is, however, not straightforward to rely on 
these studies to discard the time lag issue; given the large 
number of parameters involved in determining how long 
a merger signature will remain visible, it is plausible that 
several late-stage mergers could have been missed. Al- 
though a merger between gas-rich galaxies can leave spec- 
tacular features for a long time, viewed from the wrong 
orientation they can be completely unnoticeable. 

While it is difficult to assess the relevance for the 
timescale issue of major mergers being overlooked, we 
can be reasonably confident that the remnant will not 
look like a disk. Spheroidal and bulge-dominated galax- 
ies are usuall y said to be formed as a result of major 
mergers (e .g.. lToomrd[l977t [B arncs fe Hernquistl 119961: 
ICox et all 120061 ). However, it has also been stated 
that disks can survive some majo r mergers, especially 
if the progenitors are gas-rich (e.g . . iBarnes fe Hernquistl 
[19961: ISpringel fe Hernquistl 120051: IHopkins et aTT feoOQ). 
nonetheless these kind of merger remnants have been ar- 
gued to not lead to a large bulge grow th and significant 
BH fueling ( |Hopkins"fe Hcrnouist 2009[). Likewise, it has 
been argued that som e gas-rich mergers can lead to the 
regro wth of the disk ([Hopkins et al.l 120091 iBundv et al.l 
I2010D . Even so, the timescales involved for such a process 
can be as much as an order of magnitude larger than the 
typical quasar lifetim e of 1- 100 Mvr (e.g.. | Porciani et al.l 
l200i IHopkins et al.l 121)0511 IShen et ai1l2007l) . 

For the significant fraction of AGN hosted by disks 
found from our classification, we could safely say that 
the mechanism responsible for triggering those AGN was 
not a past major merger, suggesting also that since z ~ 1 
alternative fueling met hods seem to play a larg er role 
than usually expected. Georga kakis et al.l J2009), from 
a sample of X-ray-selected AGN, compared the lumi- 
nosity function of their disk-hosted AGN against the 
analytic model of the X-ray AG N luminosity function 
for a s tochastic accretion mode by IHopkins fe Hernquistl 
(2006). They showed that the model can reproduce the 
observations, but at the same time the overall number 
density of the observed disks was underpredicted, espe- 
cially at high X-ray luminosities. On our sample of 140 
AGN, 18 sources have Lx > 10 44 erg s _1 , from which 10 
of their host galaxies were classified as disk dominated 
with an agreement >80%. This suggests that alterna- 
tive BH fueling methods (i.e., those that do not destroy 
the disk) are not only more common on the overall AGN 
population at z < 1, but also much more efficient than 
the existing models predict. 

5.1.2. No veiled X-ray activity in merging galaxies 

The aforementioned models leave the possibility that 
we could be missing an important fraction of AGN due to 
gas and dust obscuration when a gas-rich major merger is 
taking place. Even though obscured AGN can still be de- 
tected through their hard X-ray emission (Hopk ins et alj 

6 For example, see [Schawinski et al. (2010) who make an exten- 
sive case of the time lag scenario. They propose an all-merger- 
driven AGN phase with a time lag of ~500 Myr for their sample 
of early-type galaxies at z ~ 0.05. Even if their result is mainly 
based on the interpretation of their data as a causal sequence of 
events (and is subject to alternative explanations), they caution 
that their particular sample only accounts for a very small fraction 
(~10%) of the overall AGN population found in the local universe. 



I2005a|) . it is possible that less luminous and highly ob- 
scured AGN lie below th e detection threshol d used to 
build the X-ray catalogs (|Treister et al.ll2004D . The X- 
ray properties of such obscured objects have been suc- 
cessfully studied in the literature b y the means of a 
stacking analysis of X-ray data (e.g.. iDaddi et alll2007l ; 
iFiore et all [20091) . If obscured AGN are being missed, 
they should be preferentially found in merging galaxies. 
Therefore, in order to test this scenario and search for 
this potentially buried X-ray activity, we stack all the 
inactive galaxies regarded as highly distorted. Eighty- 
seven inactive galaxies fulfill our simultaneous criteria 
of being individually classified as either Dist-1 or Dist-2 
with an agreement of >75%, and classified as Dist-2 with 
an agreement > 65%. 

For this analysis, we take advantage of the higher sen- 
sitivity of the Ch andra observatio ns of the COSMOS 
field (C-COSMOS: lElvis et al.1l2T)09l ). compared with the 
XMM-Newton data. Even though Chandra covered only 
half of the field (~0.9 deg 2 ), it has a flux limit three times 
below the XMM-Newton sensitivity, which makes the 
tradeoff in smaller coverage absolutely justifiable, con- 
sidering that we want to detect possible X-ray sources 
below the XMM-COSMOS catalog sensitivity threshold. 

For the stacking of the X-ray data, we used the 
CSTACK tooQ developed by one of the authors (T.M.), 
which includes a detailed bootstrapping error analysis 
through 500 realizations. Because the stacking is made 
from multiple observations, we consider the counts within 
a radius varying according to the off-axis angle, corre- 
sponding to 90% of the encircled counts. We stacked 
the 45 objects that lie within the C-COSMOS area, af- 
ter excluding 1 object that was close to an X-ray source. 
We found an excess of soft 0.5-2 keV and hard 2-8 keV 
count rates from the source region at modest levels of 
2.2cr and 2.4<r, respectively. Figure [8] shows the results 
of the stacking in the two energy bands, with the aver- 
age radii of 3'.'4 and 3'.'7 for comparison, within which no 
source is noticeable above the background level. 

The lack of any obvious source after the stacking sug- 
gests that this moderate excess could be in part due to 
the expected emission from star-forming galaxies, and 
also from extended source emission (e.g., from a galaxy 
group). The possibility that a few sources dominate the 
overall count rate is unlikely since (1) the shape of the 
count rate distribution is that of a unimodal Gaussian 
and (2) no outliers are present. Therefore, it is doubt- 
ful that we are missing a significant fraction of accreting 
BHs hidden within the population of inactive galaxies 
undergoing interactions. 

5.1.3. No Enhanced Soft X-ray Absorption in Merging AGN 
Host Galaxies 

As mentioned before, AGN obscuration due to the sur- 
rounding gas and dust during a major merger would af- 
fect mainly the soft X-ray energy band, while the hard 
band would remain unobscured. If we observe an AGN 
hosted by a merging galaxy, and this interaction was re- 
sponsible for the BH activity, we would expect to observe 
a hard X-ray spectrum from this source. To trace the 
obscuration level of our interacting AGN host galaxies, 
we compute their X-ray hardness ratio (HR) . The HR is 

7 http://cstack.ucsd.edu 
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Figure 8. Stacked Chandra images of 45 inactive galaxies likely 
to be undergoing a major interaction, on the soft 0.5-2 keV (left) 
and hard 2-8 keV (right) energy bands, showing the average radii 
of the stacked sources as white circles. The cutouts are 12" X 12". 

defined as 

HR={H - S)/{H + S), (4) 

where H and S stand for the hard (2-10 keV) and soft 
(0.5-2 keV) counts, respectively. At our redshift range, 
it is still safe to say that the HR values lower than -0.2 
corre spond to an unabsorbed, soft spectrum ([Hasingerl 
[HI). 

From our visual analysis, we have 13 AGN host galax- 
ies regarded as highly distorted with high agreement ac- 
cording to the criteria used before. By computing the HR 
for these objects, we find that, contrary to what models 
predict, all of these particular sources present soft X-ray 
spectra. All of them have HR values <-0.2, with a mean 
of -0.53, which shows a low attenuation in the soft band. 

It has been argued, however, that the HR diagnostic is 
rather crude in terms of predicting obscuration, and in- 
deed, bright Compton-thick AGN can feature soft X-ra y 
spectra due to photoionized gas (jLevenson et al.l f2006). 
Even so, this is only valid when the AGN is not observed 
directly, and we can easily establish that at least for the 
type-1 subsample this would not the case, and that we 
are certainly looking at active, accreting BHs. Looking 
only at the seven type-1 objects from these likely merg- 
ing galaxies, we find that the average HR is -0.56 which 
indeed suggests a low level of obscuration. 

One possible interpretation is that these interactions 
are not related to the observed AGN episode and that are 
instead only chance encounters. Dissipationless or gas- 
poor mergers could account for the lack of obscuration, 
but then it is unlikely that any strong merging signatures 
and substantial accretion onto the central BH would take 
place directly due to these kind of events. iPierce et al.l 
(2007) found the same result for X-ray-selected AGN 
hosted by interacting galaxies, suggesting that the ob- 
served interactions were not responsible for the fueling 
of those accreting BHs. 

From another perspective, however, the models men- 
tioned earlier are limited by the proposed picture 
of the merger-ULIRG-feedback-quasar ti meline (e.g., 
iSanders et allll988allb1 : iHopkins et al.ll2008f ). which is al- 
ready regarded as oversimplified. The AGN phase is said 
to happen after coalescence, but observations of large 
samples of ULIRGs, all of them undergoing interactions, 
have found a significant scatter in the trends of AGN 
contribution, accretion rate, and d ust obscuration with 
merging state (jVeilleux et al.1120091 ). Some of these have 



even been found to be dominated by the AGN in pre- 
merging state. Chaotic behavior during a merger event 
can lead to various unpredictable episodes of starburst 
and nuclear activity. Such episodic behavior can start 
much earlier than the final coalescence and can be re- 
sponsible for different periods of gas inflows, obscuration, 
and visibility, therefore explaining an already unobscured 
merger-induced AGN still early during the interaction, 
as traced by the soft spectra observed in our interact- 
ing active galaxies. This conclusion, at the same time, 
contradicts the alternative time lag scenario. 

5.2. Major Merging: Not the Most Relevant Mechanism 

Our analysis has demonstrated that the scenario in 
which mergers are responsible for triggering AGN after a 
significant time lag is unlikely. The high fraction of disks, 
the lack of a hidden significant AGN signal in merging 
inactive galaxies, and the missing soft X-ray obscura- 
tion of interacting AGN hosts all appear to rule out this 
model as a possible explanation of our results. The ab- 
sence of any further evidence in support of this scenario 
leads us to the only remaining possible interpretation of 
our results: active galaxies are involved in major mergers 
no more frequently than inactive galaxies, and mergers 
have not played a leading role in AGN triggering for the 
last 7.5 Gyr. Our results agree with the few recent stud- 
ies that have used a con t rol sample (iDunlop et al.l [20031: 
i Grogin et al.ll2003l [20051: IPierce et alJl2007t iGa bor ■ et al l 
l2009UReichard et alJl2009l : iTal et al.ll2009l and also with 
recent results from the E-CDFS by Bohm et al. 2010, 
submitted), in the sense that the morphologies of the 
AGN host galaxies are not unusual and do not show 
a preference for merging systems. Of the studies men- 
tioned earlier which supported a merger-AGN connec- 
tion, many only provided circumstantial evidence for 
such a link, without any control sample comparisons. 

The lack of enhancement on merging signatures for 
AGN hosts with respect to the background level indi- 
cates that there is no causal connection between merging 
and AGN triggering up to z ~ 1 and M* ~ 1O 1:L7 M , 
the galaxies dominating BH growth at these redshifts. 
It is still a plausible scenario that major mergers could 
be responsible for some of the brightest quasars; we do 
not intend to neglect this possibility, but in the con- 
text of a clean, large X-ray selected population of AGN, 
it is certainly not the most relevant mechanism. The 
large fraction of AGN hosted by disk-dominated galaxies 
shows that alternative mechanisms, i.e., stochastic pro- 
cesses and minor mergers dominate, for this sample of 
objects. 

The merger-starburst connection has also been widely 
studied in the same perspective. Both mechanisms share 
the need for enough cold gas to be brought to the cen- 
tral regions of the galaxy, so it is worth mentioning 
analogous conclusions from the recent literature: (1) in- 
deed, major mergers can t rigger strong s tarbursts (e.g., 
iMihos fc Hernquistl [l99l iSpringell [2000). but (2 ) not 
always, as se en in models ([Di M attco et al. 2007]) and 
observations (jBergvall et al.l I2003]). and (3) its overall 
contribution is re latively modest (|Di Matteo et al.l[2008t 
Uogee et"aT1l2009D . with no more than 10% of star forma- 
tion in mas sive galaxies being tr iggered by major mergers 
at z ~ 0.6 (jRobaina et al ll2009D. 

Different studies (e.g. JBallantvne et al.ll2006l : lHasingerl 
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[2008tlLi et aLl [2010) have converged on proposing the fol- 
lowing scenario: the major merger-driven evolution dom- 
inates early in the universe, producing the bulk of the 
brightest quasars at z = 2 — 3. Around z ~ 1 however, a 
different evolutionary mechanism takes over, with secular 
processes becoming the main triggers for the BH activ- 
ity and growth. While our analysis cannot be performed 
at higher redshifts with the current observational data 
set, our results appear to fit this picture. Nevertheless, 
the overall relevance of major merging, even in the early 
universe, has yet to be determined. Other recent studies 
suggest that secular processes play a much larger role: 
observations of massive star-forming galaxies at z ~ 2 
have shown that their buildup has been dominated by 
cold rapid accretion and secular processes (|Genzel et al.l 
2008), without the need of major mergers. It has been 
stated on the basis of dark matter simulations that the 
likely number of major mergers is insufficient to account 
for the transformation of star-for ming turbulent dis ks at 
z = 2 into ellipticals at z = (|Genel et al.l 120081) . A 
broad er view of the accre tion history of dark matter ha- 
los bv iGenel etafl (|2010t l. quantified that ~60% of the 
dark matter in a given halo is contributed by mergers, 
with only ~20% being major mergers. Instead, the rest 
(~40%) of the dark matter would be accreted smoothly. 
This also agrees with recent work using smooth par- 
ticle hydrodynamic simulations, stating that galaxies 
have acquired most of t heir baryonic mass throu gh the 
cold mode of accretion (|Kere§ et al.1 12005L 12009( 1. Fur- 
thermore, merger-free models have shown that isolated 
galaxies can reproduce the quasar duty cycles between 
z = 1 and 3 and fee d their BHs with the re cycled gas 
from evolving stars (jCiotti fc Ostrikerl 120071) and even 
repro duce the observed scaling relations (jLusso fc Ciottil 
2010). Overall, these studies have shown that secular 
evolution can be highly relevant, also at the redshifts at 
which the peak of quasar activity occurs. 

6. CONCLUSIONS 

In this work, we performed a consistent visual analy- 
sis on the ifST-based morphologies of a sample of 140 
X-ray-selected AGN host galaxies over z ~ 0.3 — 1.0 and 
M* < 10 1 M Q , and compared them with a matched 
control sample of inactive galaxies under the same con- 
ditions. Our goal was to search for the presence of any 
significant connection between major merging and BH 
fueling as suggested by models and observational tests. 
In summary: 

1. From our visual analysis, ~85% of our AGN host 
galaxies show no strong distortions on their mor- 
phologies. Comparison with the control sample 
shows that the distortion fractions are equal within 
the 68% central confidence level. Given our sample 
size, we can state that at a 95% confidence level 
the highly distorted fraction of AGN hosts is less 



than 1.9 times that of the inactive galaxies. Merg- 
ers and interactions involving AGN hosts are not 
dominant, and occur no more frequently than for 
inactive galaxies. 

2. Over 55% of the AGN from our sample are hosted 
by disk-dominated galaxies, implying a triggering 
mechanism that would not destroy the disk, i.e., 
not a major merger. This also indicates that it 
is unlikely that we could be missing major merg- 
ers due to strong distortions having already been 
washed out over a large time lag prior to the ig- 
nition of the AGN. The presence of an impor- 
tant fraction of disk-dominated hosts on the AGN 
brighter than Lx > 10 44 erg s _1 suggests that sec- 
ular fueling mechanisms can be highly efficient as 
well. 

3. Through a detailed stacking analysis of the X-ray 
data of our inactive galaxies undergoing mergers, 
we did not find an underlying X-ray signal indi- 
cating the presence of a substantial population of 
obscured AGN. 

4. Looking at the hardness of the X-ray emission of 
our AGN hosts that are clearly undergoing an inter- 
action, we found soft X-ray spectra in all of them, 
contradicting the expected obscuration in this band 
predicted by models. This can be either because 
the observed interactions are not responsible for the 
BH fueling or the unpredictable output of a merger 
event allows many accretion phases as well as an 
unobscured AGN, even during such early stages. 

Our work explicitly suggests that, at least for the last 
7.5 Gyr, major merging has not been the most relevant 
mechanism in the triggering of typical AGN, and that the 
bulk of the BH accretion occurs through internal secu- 
lar processes and minor interactions. The alternative 
interpretation of a time lag between merger trigger and 
AGN onset is unlikely due to the zero enhancement of the 
distortion fraction, the high incidence of disks, and the 
absence of a significant X-ray signal in merging inactive 
galaxies as a potential buried AGN population. 
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APPENDIX 

A. AGN-HOST GALAXY DECOMPOSITION 

The light distribution of the type-1 AGN is clearly dominated by the bright active nucleus, and because we want 
to analyze the morphologies of their host galaxies, accurate removal of the nuclear source is of vital importance. 
This is done through a rigorous two-dimensional parametric fitting with GALFIT, with which we reduce each system 
down to a two-component model: a PSF to represent the AGN and a Sersic light profile accounting for the host 
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Figure 9. The difference in the observed magnitudes (If814,w) of the comparison galaxies before (in) and after (out) the point source 
addition/subtraction. The left-hand panel plots this difference against the initial magnitude, and the right-hand panel against the host to 
nucleus flux ratio, H/N. The 1<t deviation away from the mean is 0.23 mag, indicated by the shaded area centered at 0.03 mag. 

galaxy. After subtraction of the modeled PSF, we are left with the host galaxy emission plus some residuals. Previous 
simulations have shown that, at our resolution and S/N, it is sufficien t with a single-component model to accoun t 
for the host galaxy rather than a more complex, multi-component one ([Sanchez et al.ll2004i ISimmons fc Urrvl[200l . 
Regarding double-nucleus sources, which could nee d a dedicated model ing, we only found 1 object in our sample. This 
is also consistent with the visual inspection from ICivano et all 1)20101 ) on C-COSMOS sources. We checked for this 
particular source (Figure [31 bottom left), for which our method removes the brightest of the two nuclei. The other 
point source is significantly fainter and does not dominate the overall galaxy brightness, hence not requiring a more 
complex decomposition. 

An appropriate initial guess of t he parameters is recomm ended to get a faster and converging model with GALFIT. 
We opt to run Source Extractor ([Bertin fc Arnouti IT996) on our cutouts to generate, in a fast and automatic way, 
rough estimates of the free parameters of the Sersic profile, such as coordinates within the image, observed magnitude, 
axis-ratio b/ a, half-light radius R e , and position angle. 

To ensure a reliable decomposition, we take particular care that the host galaxy is modeled with the least possible 
unnecessary flux transfer between PSF and the Sersic profile, which would result in either an over- or undersubtraction. 
We perform several GALFIT runs on each obj ect with 3 different ch oices for the Sersic index n: we fix it to a n = 1 
exponential profile (jFreemanlll970|) . to a n — 4 lde Vaucouleursj (jl948[) profile, and we also leave it a as free parameter 
for GALFIT to decide. To choose the right model, we need our best fit to be reasonable in terms of the resulting 
parameters. We require our host galaxy model: (1) not to be too concentrated or too shallow, meaning a half-light 
radius between 2.5 pixels < R e < 100 pixels, (2) not to diverge to extreme elongations, therefore to have b/a > 0.5, 
and (3) to have its Sersic index within 0.5 < n < 8, for the free n case. We interpret that if the values run away 
from these boundaries, GALFIT did not manage to model the underlying galaxy but instead could be accounting for 
uncertainties in the PSF. 

The model with the least \ 2 is chosen between those that comply with the above criteria. If the model with the free 
Sersic index is chosen among the 3 cases, we reassign the index and model it as follows: (1) if n < 2 then n = 1, (2) if 
2 < n < 3 then n = 2.5, and (3) if n > 3 then n — 4. 

A key aspect of the AGN-host galaxy decomposition is the choice of an accurate PSF, both for modeling the AGN 
itself and for deconvolving the host galaxy light distribution. Even though the space-based HST provides extremely 
precise PSFs, instrumental effects are still important. The position of the target within the detector and the temporal 
variability along different orbits can lead to discrepancies between the PSFs from the observations and the ones used 
for the analysis. This yields systematic errors in the image decomposition which can be critical for very bright AGN. 
The COSMOS survey provides us with the opportunity to minimize these spatial and temporal effects by using stellar 
PSFs from stars observed under the same conditions than our targets. For each object, we build spe cific PSFs by 
avera ging the neares t ~30 stars in the same manner to other similar studies with large HST coverage (jJahnke et al.l 
2004; Sa nchez et al"][200l . The rms error from the creation of each PSF is propagated to the intrinsic variance of the 
AGN; the uncertainty of the object being fitted is required by GALFIT in order to converge to a minimum normalized 

x 2 - 

B. CREATfNG MOCK AGN HOSTS 

For our subsample of type-1 AGN, we build a special comparison sample of simulated AGN hosts by adding stars as 
fake nuclei to our inactive galaxies. To remain true to the characteristic blue colors of the AGN, we perform an initial 
selection of stars from the COSMOS ACS archive by placing color cuts in (B - V) < 0.75 and (V - R) < 0.95. For 
each of the control galaxies, we look for stars that match the contrast level between the fluxes of the host and nucleus 
(H/N) of the corresponding AGN. With a matching star found, we simply add it over the centroid of the galaxy. 

We then apply the same point source removal procedure as for the original type-1 sources. PSFs are created exactly 
as before, and the light contributions of the star and the underlying galaxy are separated using GALFIT. With the 
exception of 3 unsuccessful fits, we are left with 727 simulated nucleus-subtracted AGN host galaxies that will serve 
as an appropriate comparison sample for our type-1 AGN hosts. 
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C. TESTING THE RELIABILITY OF THE IMAGE DECOMPOSITION 



The creation of a sample of simulated nucleus-subtracted hosts from a starting point of real galaxies and stars gives 
us the opportunity to check the impact of our point source removal technique, and to see whether this technique is 
biased. How significant are the residuals? We have performed photometry on the control galaxies before and after the 
addition/subtraction of the fake nucleus. If a large magnitude offset were to be found, we would have had to consider 
reselecting our control sample, because we would inevitably be comparing active and inactive galaxies with different 
observed magnitudes. We find that, on average, the galaxies are fainter by 0.03 mag after the subtraction, with a la 
deviation of 0.23 mag. Figure [9] shows the difference between the initial and recovered magnitudes for the hosts as a 
function of the initial magnitudes and H/N ratio for our control galaxies. There is no obvious correlation between the 
offset and the initial magnitudes of the galaxies, but as expected the recovered values tend to be less exact for more 
compact galaxies and brighter active nuclei. 

These results show that this technique is trustworthy, and the offset found can be considered negligible and does not 
affect our choice of a comparison sample. 
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